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ABSTRACT 

The far-infrared and submillimeter portions of the electromagnetic spectrum provide a unique view of the astro- 
physical processes present in the early universe. Our ability to fully explore this rich spectral region has been 
limited, however, by the size and cost of the cryogenic spectrometers required to carry out such measurements. 
Micro-Spec (/i-Spec) is a high-sensitivity, direct-detection spectrometer concept working in the 450-1000 /im 
wavelength range which will enable a wide range of flight missions that would otherwise be challenging due to 
the large size of current instruments with the required spectral resolution and sensitivity. The spectrometer de- 
sign utilizes two internal antenna arrays, one for transmitting and one for receiving, superconducting microstrip 
transmission lines for power division and phase delay, and an array of microwave kinetic inductance detectors 
(MKIDs) to achieve these goals. The instrument will be integrated on a ^ 10 cm^ silicon chip and can therefore 
become an important capability under the low background conditions accessible via space and high-altitude 
borne platforms. In this paper, an optical design methodology for //-Spec is presented, with particular attention 
given to its two-dimensional diffractive region, where the light of different wavelengths is focused on the different 
detectors. The method is based on the maximization of the instrument resolving power and minimization of the 
RMS phase error on the instrument focal plane. This two-step optimization can generate geometrical configu- 
rations given specific requirements on spectrometer size, operating spectral range and performance. Two point 
designs with resolving power of 260 and 520 and an RMS phase error less than ^ 0.004 radians were developed 
for initial demonstration and will be the basis of future instruments with resolving power up to about 1200. 
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1. INTRODUCTION 

Far- infrared (IR) and submillimeter (15 fim-1 mm) spectroscopy provides a powerful tool to probe a wide range 
of environments in the universe. In the past thirty years, discoveries made by several space-based observatories 
have provided unique insights into physical processes leading to the evolution of the universe and its contents. 
This information is encoded in a wide range of molecular lines and fine structure lines. Observations of such 
spectral lines enable the exploration of galaxies at high redshifts. The fine structure lines of abundant elements 
(C, N, and O) allow tracing the obscured star formation and AGN activity into the high-redshift universe. One 
can measure galaxy redshifts and determine their elemental abundances and physical conditions out to redshifts 
of 2 ; > 5. In spite of this, a number of questions remain unanswered regarding the very early steps of the universe 
as well as galactic, stellar and planetary formation. The ability to fully explore this rich spectral region has been 
limited by the size and cost of the cryogenic spectrometers required to carry out these measurements. The work 
proposed here specifically addresses the need for integrated spectrometers and ultra-low-noise, far-IR, direct 
detectors, whose specific requirements are shown in Table and compared against the current state of the art. 

To contribute to realizing the goals outlined in Table 1 and fill in the gaps in the current state of the art, 
a high-performance integrated spectrometer module, Micro-Spec (//-Spec), operating in the 450-1000 //m (300- 
650 GHz) range is proposed. With //-Spec (Fig. 1), the incoming radiation collected by the telescope is coupled 
to the spectrometer via a broadband slot antenna and directed to a series of power splitters and a delay network 

Further author information: (Send correspondence to G.C.) 

G.C.: E-mail: Giuseppe.Cataldo@nasa.gov, Telephone: -hi 301 286 7497 


Table 1: Summary of far-IR spectrometer and detector array requirements and comparison with current state of 
the art. 


Technology 

Metric 

Symbol 

State of the art 

Requirements 


Wavelength range 

A 

250 — 700 jam 

220 - 2000 


Sensitivity 

NEP 

10-19 W/VHz 

3 X 10-91 W/VHz 


Resolution 

n 

> 100 

> 1200 

Far-IR detector arrays 

Detective Quantum Efficiency 
Format 

DQE 

~ 15% 

30 X 30, 

Component test 

> 90% 

1500 spectral 
elements 


Temperature 

T 

< 1 K 

- 0.05 K 


Time constant 

T 

100 ms 

< 10 ms 


made of superconducting microstrip transmission lines. The delay network creates a retardation across the input 
to a planar waveguide multimode region, which has two internal antenna arrays, one for transmitting and one 
for receiving the radiation as a function of wavelength. Absorber structures lining the multimode region control 
the power emitted into large angles or reflected from the receiver antenna array. An array of feed structures is 
employed to couple the radiation to the multimode region and concentrates the power along the focal surface with 
different wavelengths at different locations. The outputs are connected to a bank of order-sorting filters which 
terminate the power in an array of microwave kinetic inductance detectors (MKIDs) for detection and readout. 
The entire instrument is integrated on a 100-mm-diameter silicon chip. This size in reduction is accomplished 
by using single-mode microstrip delay lines, which can compactly be folded on the silicon wafer by a factor of 
the medium’s effective refractive index. 

Filter MKID 

i+1 
i 

i-1 


Figure 1: Layout of the /i-Spec module. The radiation is coupled into the instrument through a broadband 
antenna and is then transmitted through a superconducting transmission line to a divider and a phase delay 
network. The spectrum enters the multimode region through an array of feeds which concentrates the power 
along the focal surface with different wavelengths at different locations. The receivers are connected to a bank 
of order-sorting filters and MKID detectors.^ 

It is worth noting that the frequency range of the implementation presented here is limited to wavelengths 
A > 250 jam by the gap frequency of currently available low-loss superconductors. These include niobium (Nb) 
and niobium-titanium nitride (NbTiN) for the transmission line structures, and molybdenum nitride (MoN) for 
the detectors. 

Finally, /i-Spec can be compared to a grating spectrometer,^ in which a plane wave is reflected from the 
grating and the phase of each partial wave scattered from the rulings is a linear function of position across the 
grating. However, it differs from similar technologies by the order of processing of the light in the spectrometer. 
For instance, in a Rowland spectrometer the required phase retardation is generated by reflection from the 
grating grooves, whereas in Z-Spec, which is an example of planar Rowland grating architecture, propagation 
occurs in parallel-plate waveguides. A last comparison can be made with one-dimensional bootlace lenses 
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found in microwave practice/^^^^ which /i-Spec builds upon for submillimeter wave applications. This paper 
will describe the design process of the multimode region and illustrate the results in terms of geometry, imaging 
quality and efficiency. 

2. MULTIMODE REGION DESIGN IN HIGHER ORDER 

A prototype version with a resolving power 7^ = 65 in first order was designed^ and built, and is currently under 
evaluation at NASA’s Goddard Space Flight Center. This design was specified by requiring the phase error to 
vanish at three specific stigmatic points on the focal plane. This can be seen as a generalization of the approach 
described in Ref. 8 with two stigmatic points. Here, we present designs for IZ = 260 (Design A) and IZ = 520 
(Design B) in higher order with zero phase error at three preselected stigmatic points employing a multi-step 
optimization process described in the following section. 

2.1 Problem formulation 

As explained in Ref. 2, the design variables are the x and y coordinates of the emitters’ centers and the 
electrical path lengths in silicon, for each feed’s electrical delay. The resolving power is defined as 

n = M-Ne, ( 1 ) 

where M is the order of the spectrometer and Ne is some power of 2, given the structure of the power divider 
network (see Fig. 1). The first step of the design consists of finding the maximum achievable resolving power, 
'TZjnax^ as a function of M and the relative emitter pitch, y = p/ X b {p = emitter pitch, Xb = wavelength associated 
with the central stigmatic point on the x axis - see Fig. 2), given specific requirements on spectrometer radius, 
R, and operating spectral range as well as certain constraints on performance. 
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Figure 2: Simplified representation of the grating geometry. On the left side three radiators can be seen, which 
point to the blaze point, B. The path from each radiator’s phase center to the first and second stigmatic points 
are also indicated by solid lines for the array’s and central reference feed. 


The formulation of this mixed integer non-linear problem is as follows: 




mdiXlZmaxiM, 77 ) = M /77 • R/^B 

( 2 ) 

subject to He{M,T]) 

< 

R (width of emitter array < radius, so that aberration is OK) 

( 3 ) 


< 

R (width of receiver array < radius, so that aberration is OK) 

( 4 ) 

Bmax ( Af, 77 ) 

> 

64 (maximum resolving power > 64) 

( 5 ) 

p(M, 77 ) 

= 

s(M, 77 ) (emitter pitch equal to receiver pitch) 

( 6 ) 

77 

> 

0 

( 7 ) 

1 < 

M 

<10, M integer 

( 8 ) 


Here we note that constraint ( 6 ) was used to simplify the positioning of the receivers on the focal plane and may 
be removed in future designs. The objective spaces as a function of M and 77 for the two cases, IZ = 260 and 
IZ = 520, are shown in Figs. 3-4, respectively. On the left (Figs. 3a-4a), it is possible to visualize the feasible 
objective space of the optimization problem described above for IZmax along with the active constraints, Eq. (4) 
(blue area) and Eq. ( 6 ) (black lines). The optimal solution to Problem (2) is thus given by the intersection 
of the contour plot representing the feasible values of IZmax and constraint ( 6 ) (lower black line). Table 2 
shows the requirements on spectrometer size and spectral range for the two cases. The minimum and maximum 
frequencies are associated with stigmatic points 2 and 1 respectively, whereas the average frequency, computed 
as their geometric mean, is associated with the blaze point, B (Eig. 2). Table 3 shows the values of the design 
variables associated with the optimal solutions as well as the values of the constraints, which are all satisfied. 

Once this problem solved, it is then possible to compute Nq = TZmax/M and round this down to a power of 
2 . According to Eq. ( 1 ), this causes the actual resolving power, calculated with this updated value of Ne, to be 
lower than IZmax- The plots in Eigs. 3b-4b show the values of the actual resolving power IZ (red line) and the 
values of M and 77 that would make such realizations possible. 

Table 2 : Requirements on spectrometer size and spectral range 


Parameter 

Symbol 

Unit 

Value (Design A) 

Value (Design B) 

Multimode region radius 

R 

cm 

1.25 

1.25 

Minimum frequency 

fmin 

GHz 

570.0 (in order M' = 

1) 

605.0 (in order M' = 4) 

Maximum frequency 

fmax 

GHz 

650.0 (in order M' = 

1) 

650.0 (in order M' = 4) 

Average frequency 

favg 

GHz 

608.7 (in order M' = 

1) 

627.1 (in order M' = 4) 

Minimum frequency 

fmin 

GHz 

142.5 (in order M' = 

4 ) 

302.5 (in order M' = 8 ) 

Maximum frequency 

fmax 

GHz 

162.5 (in order M' = 

4 ) 

325.0 (in order M' = 8 ) 

Average frequency 

favg 

GHz 

152.2 (in order M' = 

4 ) 

313.5 (in order M' = 8 ) 


Table 3: Optimal solutions of Problem ( 2 ) for designs A and B 


Variable 

Symbol 

Unit 

Value (Design A) 

Value (Design B) 

Spectrometer order 

M 

- 

4 

8 

Emitter diameter 

77 

- 

1.1634 

1.1809 

Maximum resolving power 

V 

f^max 

- 

297.83 

604.60 

Emitter array width 

He 

cm 

1.0912 

1.0751 

Receiver array width 

Hr 

cm 

1.1297 

1.2154 

Emitter /receiver pitch 

p = s 

cm 

0.0165 

0.0163 

Number of emitters 

Ne 

- 

65 

65 

Number of receivers 

Nr 

- 

87 

91 

Resolving power 

■R 

- 

260 

520 



(a) The contour plot represents all feasible values of (b) The contour plot represents the feasible values of 

IZmax and shows a solution exists for M = 4 that sat- IZ for all orders given as powers of 2. The IZ — 260 red 

isfies all constraints. curve identifies the optimal solution. 

Figure 3: Objective spaces of Problem (2) for Design A. The active constraints are Eq. (4) (blue area) and 
Eq. (6) (black lines). 
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(a) The contour plot represents all feasible values of 
IZmax and shows a solution exists for M = 8 that sat- 
isfies all constraints. 

Eigure 4: Objective spaces of Problem (2) for Design B. 
(black lines). 


M 

(b) The contour plot represents the feasible values of 
IZ for each order given as powers of 2. The IZ = 520 
red curve identifies the optimal solution. 

The active constraints are Eq. (4) (blue area) and Eq. (6) 


The second step consists in minimizing the overall RMS phase error, (Prms^ on the focal plane in order to 
determine the optimal solution in terms of the above-mentioned design variables. 


min 



^RMS d6> = 



[iPij{xi,yi,Ti,Oj) - {(fjOjW 
Ne 


d(9 


{ iPii{xi,yi,Rf,Oi) = 0 

(fi2{xi,yi,Rl,02) =0, i = 1, A^e 

Vii •) ^ 0 


(9) 

(10) 


Here, ipij is the relative phase of each transmitter, {(f{0j)) = 0 is the relative phase of the central transmitter (this 
is zero by construction as the central radiator is used as a reference) and Oj represents the angle corresponding to 
each of the points in which the focal plane was discretized, ^i, O 2 and Or are associated with the three stigmatic 
points as depicted in Fig. 2. 


2.2 Optimization results 

The solutions to the system of equations (9)- (10) found for the two cases discussed in this paper can be seen in 
Fig. 5. The emitters’ positions are indicated in red and present several characteristics which are worth discussing. 
First, they do not lie on the grating circle but on a curve that is tilted leftwards and intersects the grating circle 
at the central emitter before ending up inside the multimode region. It was verified that the shape of this curve 
only approximates a circle with a radius ^ 2.2R, is not symmetric and is caused by the imposition of zero 
RMS phase error on the blaze point, B (last constraint in Eqs. (10)). It was verified that if such a constraint 
is removed, the emitters lie exactly on the grating circle, as shown in the literature.^’ ^ In our previous work,^ 
this tilting effect was explicitly avoided by constraining all the emitters to lie within Ai/8 from the 2R circle. 
Figure 5b shows also that it is less accentuated for the configuration in order M = 8. 




(a) M = 4 and 7^ = 260 (b) M = 8 and = 520 

Figure 5: Multimode region design for the two cases discussed above. 

The RMS phase error is shown in Fig. 6 for the two cases presented in this paper and compared with that 
computed for the optimized IZ = 6b version described in Ref. 2. The new designs represent an improvement of 
no less than ^ 64% over the nominal spectral range indicated in Table 2. While the RMS phase error vanishes at 
the three imposed stigmatic points, it is worth noting the presence of a fourth stigmatic point for both designs at 
^ 130° and ^ 155°. It can be shown that, in general, the RMS phase error as defined in Eq. (9), is a fourth-order 
function of the focal plane angle, 6>, and thus admits a maximum of 4 zeros. Here, equality constraints such 
as those in Eqs. (10) were only imposed on three points, leaving the fourth free because beyond the spectral 


range of interest. As a consequence, this does not affect the overall quality of the instrument for the specified 
requirements. Exploration of a configuration with 4 stigmatic points spanning a desired frequency range and 
such to minimize the RMS phase error as defined in Problem (9) is left for future effort. 



Figure 6 : RMS phase error for the two designs discussed in this paper and compared to the optimized 7^ = 65 
version presented in Ref. 2 . 

2.3 Power efficiency 

The power coupling in the new design configurations was computed with the model described in Ref. 2 , Sec. 4. 
Figure 7 illustrates the values of the normalized power for the IZ = 260 run in first, second and fourth order, 
as well as for the IZ = 520 run through eighth order. When comparing the coupling efficiency of these designs 
in first order with that described in Ref. 2 , it can be seen that the TZ = 260 version outperforms the others, 
though still inferior to the optimized 7Z = 65 version. This low efficiency is the result of higher-order diffraction 
peaks that show up in the multimode region and are caused by 77 being greater than 1 / 2 , as required to avoid 
diffraction.^^ 


3. CONCLUSIONS 

A design methodology has been developed for high-resolution configurations of the /i-Spec multimode region. 
The design procedure first maximizes the resolving power subject to constraints on geometry, operating frequency 
range and performance, thereby determining the order of the spectrometer. This then allows the RMS phase 
error on the instrument focal plane to be minimized and vanish at three preselected stigmatic points associated 
with three specific wavelengths. The two designs explored here achieve resolutions of 260 and 520 in fourth 
and eighth order respectively and present a coupling efficiency approaching unity. The coupling efficiency was 
observed to decrease monotonically with the order. The RMS phase error between the three stigmatic points is 
less than 0.005 radians for both designs and the presence of a fourth stigmatic point was observed beyond the 
angular range in use. This feature can be used in future work to increase the number of spectrometer channels 
and resolving power. 
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Figure 7: Comparison of the efficiency of the two designs presented in this paper and the ones from Ref. 2. 
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